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SUMMARY 


PROBLEM 

Analyze  the  vertical  directionality  of  the  ambient  noise  in  the  Ionian  Basin  of  the 
Mediterranean  Sea.  Record  and  analyze  data  for  directional  content  using  adaptive  and 
conventional  beam  forming  techniques. 

RESULTS 

It  was  found  that  the  amount  of  uncorrelated  (circuit)  noise  has  a  strong  influence 
on  the  measured  directionality  of  the  noise  field.  It  was  also  found  that  adaptive  beamform¬ 
ing  improves  the  resolving  power  of  the  array  at  any  frequency  over  conventional  beam¬ 
forming.  but  more  especially  at  low  frequencies. 

RECOMMENDATIONS 

1 .  Array-independent  vertical  directionality  of  ambient  noise  should  be  attempted 
using  adaptive  beamforming  techniques  as  an  analytical  tool.  Minimal  ambiguity  in  direc¬ 
tionality  conclusions  will  result  from  careful  attention  to  uncorrelatcd  (circuit)  noise,  sus¬ 
pension  noise,  and  verticality  of  the  array. 

2.  The  resolving  power  of  the  array  should  be  tested  in  advance  using  propagation 
theory  and  assumed  shipping  and  sea  state  distributions. 

3.  Only  then  can  measurement  results  help  in  our  understanding  of  the  acoustic 
environment  and  allow  us  to  establish  system  performance  modeling  as  a  basis  for  future 
system  decisions  with  adequate  confidence. 
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INTRODUCTION 

(U)  In  January,  1970,  under  contract  N00019-70-C-0468,  extensive  data  were  taken 
in  the  Mediterranean  Sea  by  Sanders  Associates,  Inc.  of  Nashua,  N.  H.  Data  were  recorded 
from  each  of  the  four  major  deep  water  basins.  The  only  one  of  these  samples  treated  in 
this  report  was  recorded  on  28  January  from  the  vertical  array  suspended  from  the  surface 
of  the  Ionian  Sea;  this  sample,  chosen  from  that  day’s  recordings,  was  recorded  during  a 
period  with  no  especially  close  merchant  traffic.  This  was  verified  by  aircraft  surveillance  as 
shown  in  figure  I. 

(U)  Figure  2  shows  the  array  as  deployed  during  the  data  acquisition.  It  consisted 
of  10  omnidirectional  hydrophones  suspended  vertically  from  a  surface  buoy  along  a  200- 
foot  cable.  These  sensors  were  recorded  as  two  separate  arrays:  one  designed  for  high 
frequencies  (300-1000  Hz)  and  the  other  for  low  frequencies  (<  300  Hz).  Data  from  the 
latter  array  were  selected  for  analysis  by  NUC.  The  sensors  of  this  array  were  regularly 
spaced  at  16.4-foot  intervals  along  the  line  configuration;  of  the  seven  sensors  recorded,  the 
elements  at  top  and  bottom  failed  to  operate  during  the  tests,  so  only  the  remaining  five 
constituted  the  array  from  which  data  were  analyzed  for  this  report. 

(U)  From  the  surface  float,  the  array  was  hardwired  to  recording  equipment  aboard 
the  U.  S.  S.  Hoist  at  the  end  of  a  2000-foot  cable  run.  This  tether  suspension  caused  numer¬ 
ous  problems,  since  the  strong  frequency  lines  from  rotating  machinery  aboard  the  Hoist 
interfered  with  fur-field  inputs.  Other  than  Hoist,  the  nearest  merchant  vessel  at  1400  on 
28  January,  the  period  taken  for  our  analysis,  was  about  20  nautical  miles  away  (see  fig¬ 
ure  1). 

(U)  The  analog  outputs  from  the  five  sensors  were  passed  through  a  30-  to  1000-Hz 
bandpass  filter  aboard  the  Hoist  just  prior  to  recording.  During  playback,  before  digitiza¬ 
tion,  antialiasing  filters  that  have  an  18  dB  per  octave  roll-off  with  the  cornering  frequency 
set  at  300  Hz  are  applied.  Digitization  was  then  set  at  1285  samples  per  second  per  channel, 
ten  significant  bits  per  sample.  With  the  data  in  the  time-domain  digital  form,  the  next  op¬ 
eration  was  1024-point  fast  Fourier  transformation  (FFT)  with  about  20  significant  bits  of 
accuracy,  using  a  computer  software  algorithm.  This  yielded  512  complex  frequency  coeffi¬ 
cients  with  approximately  1.25-Hz  resolution. 

(U)  Owing  to  the  regular  sensor  spacing  along  the  65.6-foot  line  array,  not  all  of 
these  frequency  bins  were  useful  for  analysis.  If  we  assume  5000  feet  per  second  as  the 
speed  of  sound  in  sea  water,  and  use  the  16.4-foot  sensor  separation  as  half  wavelength  for 
the  upper  frequency  of  interest,  we  find  that  the  array  performance  will  gradually  degrade 
until  at  1 52.4  Hz,  directional  ambiguity  will  prevent  adequate  resolution  of  the  directional 
character  of  the  noise  field.  Hence,  we  lack  confidence  in  any  conclusions  about  results 
above  1 52.4  Hz.  After  FFT,  the  complex  frequency  coefficients  were  multiplied  to  yield 
spectral  crosspower  matrices  by  adding  successive  sensor  pair  multiples  until  an  effective  8.5 
minutes  of  time  integration  was  achieved. 

(U)  These  matrices  were  “interrogated”  for  the  directional  character  of  the  field  by 
pre-  and  post-multiplication  by  a  directional  pointing  vector  which  describes  the  expected 
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«')  figure  1.  Surface  traffic  plot . 

wavefront  of  acoustical  energy  from  a  single  direction  and  at  a  single  frequency.  The  result¬ 
ant  integrated  beam  power  outputs,  by  repeated  interrogation  over  all  directions  and  all 
frequencies,  can  be  plotted  as  shown  in  figures  3  (adaptive)  and  4  (conventional)  to  show  an 
overall  measure  of  the  actual  data  noise-field  directionality.  For  single  frequency  considera¬ 
tion,  however,  the  plots  (figures  A 1  -  A 12)  from  actual  data  can  be  better  compared  with 
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(U)  Figure  2.  Vertical  line  array. 

a  priori  assumptions  about  the  vertical  directionality  (see  figures  B1  -  B72).  All  the  above 
computations  were  based  upon  the  same  well-established  theory  found  in  reference  3.  and 
the  reader  is  referred  there  for  more  detailed  information  about  the  beamforming  tech¬ 
niques  employed  for  this  analysis.  For  further  details  regarding  the  array  and  data  collec¬ 
tion  procedures  see  the  final  Sanders,  Inc.  report  (reference  4). 
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(U)  Figure  3.  Adaptive  FRD/Ii. 
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(U)  Figure  4.  Convention*!  FRD/E. 
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RESULTS 

(U)  Prior  to  the  analysis  of  actual  array  data,  we  tested  the  ability  of  the  array  to 
resolve  vertically  directional  simulated  noise.  To  do  so,  we  chose  the  upper  solid  curves  in 
figures  B1  through  B72  as  the  simulated  directionality.  We  then  employed  two  bcamform- 
ing  schemes  to  resolve  this  known  directional  character:  adaptive  (or  optimal)  combination 
of  the  five  sensors,  and  conventional  (or  time-shift-and-sum).  Results  of  the  application  of 
these  two  beam  former  outputs  are  shown  in  figures  B 1  through  B72  as  the  middle  and  lower 
curves,  respectively.  Note  that  at  the  lower  frequencies  (<  64  Hz)  the  resolving  power 
decreases  sharply,  as  indicated  by  the  lower  (conventional  beamforming)  curve,  which  is  no 
doubt  due  to  the  array’s  short  (65.4-foot)  length.  Note  also  that  at  frequencies  above  the 
half-wavelength  frequency  (152.4  Hz)  degradation  can  be  seen  (especially  in  the  convention¬ 
al  response)  as  a  waviness  in  the  baseline  in  directions  away  from  the  directionality  peak, 
which  becomes  progressively  worse  as  frequency  is  increased.  In  this  upper  frequency  band 
the  array  is  sensor-limited  and  the  regular  spacing  between  sensors  leads  to  “grating  lobe” 
behavior  or  directional  ambiguity  at  higher  frequencies.  Adaptive  beamforming  improves 
the  resolving  power  of  the  array  at  any  frequency,  but  is  especially  effective  at  low  frequen¬ 
cies.  Figures  5-10,  A1-A12,  and  B1-B72  present  the  data  as  “bearing  response”,  plotting 
relative  bearing  (0°  straight  up,  90°  horizontal,  and  180°  straight  down)  and  relative  power 
out  in  dB  along  the  ordinate.  No  attempt  is  made  to  indicate  absolute  power  levels,  and  the 
curves  are  arbitrarily  clamped  to  0  dB  at  180°.  This  type  of  output  is  intended  for  direc¬ 
tionality  conclusions  only:  other  outputs  (array  gain  versus  frequency)  are  more  useful  for 
array  S/N  improvement  conclusions. 

(U)  Actual  data  results  to  be  compared  with  the  above  a  priori,  simulated  data 
results  are  found  as  figures  A 1  through  A 12.  Note  that  they  all  show  broader  directional 
response  than  any  of  the  simulated  data  results  (figures  B1  through  B72).  Note  also  that  the 
beam  power  response  humps  in  the  actual  data  (figures  A 1  through  A 12)  suggest  that  the 
actual  data  field  is  best  approximated  by  either  noise  field  1  (figures  B1  through  B 12)  or  2 
(figures  B13  through  B24),  depending  upon  the  frequency  being  considered.  However,  the 
separation  between  the  adaptive  and  conventional  beam  power  output  response  curves  is  less 
for  the  actual  data  than  for  the  corresponding  simulated  data  response  curves.  This  observa¬ 
tion  caused  us  to  examine  the  conditions  from  which  this  result  would  follow.  We  conclud¬ 
ed  that  this  was  probably  due  to  high  levels  of  uncorrclated  (circuit)  noise,  and  we  present 
the  following  sample  computational  results  to  document  this  conclusion  and  to  arrive  at  a 
reasonable  estimate  of  the  vertical  angle  arrival  structure  of  the  actual  noise  field. 

Directional-Response  Computations  at  97  Hz 

(U)  In  figure  5  we  show  a  family  of  four  simulated  directional  response  curves  (the 
upper  four)  so  constructed  that  they  could  account  for  the  actual  data  directional  response 
at  97  Hz  (which  is  presented  in  figure  6  and  shown  superimposed  on  the  simulated  direction¬ 
al  response  curves  in  figure  7)  under  different  assumptions  about  the  fraction  of  uncorrclat¬ 
ed  noi;;e  on  each  sensor  channel.  In  figure  5,  the  top  curve  in  the  family  portrays  the  actual 
vertical  character  of  noise  in  the  water,  which  could  account  for  the  response  shown  in  the 
two  lowci  curves  if  0. 1  or  -  1 0  dB  of  the  single  sensor  voltage-squared  input  had  been  con¬ 
taminated  with  uncorrclated  (from  channel  to  channel)  noise  such  as  that  from  thermal  or 
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circuit-noise  sources.  The  next-lower  curve  of  the  family  is  the  constructed  response  (not 
unique)  appropriate  for  9-percent  contamination,  the  next  for  5-percent  contamination,  and 
the  lowest  family  member  for  no  contamination.  Which  curve  of  this  family  best  represents 
the  true  directional  character  of  the  vertical  noise  field  at  the  time  and  location  of  the  de¬ 
ployed  Sanders  array  cannot  be  conclusively  established,  which  underscores  the  need  for 
high  quality  array  recording  (low  circuit  noise)  if  directionality  conclusions  are  to  be  drawn 
from  such  measurements. 

(U)  A  check  with  Sanders,  Inc.,  revealed  that  the  spectral  circuit  noise  level  had 
been  measured  at  only  about  10  dU  below  the  single  sensor  power  at  the  recording  time,  so 
we  continued  our  analytical  investigation.  We  took  the  actual  data  spectral  cross  power 
matrix  at  97  11/  and  subtracted  from  the  main  diagonal  terms  a  scalar  amount  equivalent  to 
10  dB  uncorrelated  noise.  Upon  interrogation,  this  matrix  yielded  the  curves  of  figure  8, 
which  are  to  be  compared  with  the  simulated  curves  of  figure  9.  f  igure  9  shows  the  simulat¬ 
ed  beam  for  liter  output  responses  assuming  a  very  small  fraction  (0.0001  or  40  dB)  of 
uncorrelated  noise.  When  we  tried  to  subtract  12-percent  uncorrelated  noise  the  resultant 
matrix  was  singular,  for  the  above  reasons,  the  upper  solid  curve  of  figure  9  is  presented  as 
the  most  reasonable  estimate  for  the  shape  for  the  actual  noise  field's  vertical  arrival  struc¬ 
ture,  bearing  in  mind  that  the  width  of  the  response  near  horizontal  may  be  artificially 
broad,  if  the  array  had  significant  tilt.  Tilt  was  not  measured.  The  curves  of  figures  8  and  9 
arc  shown  superimposed  in  figure  10  to  facilitate  comparison. 

Array  Gain  vs  Frequency  (AGVF)  Curves 

(C)  If  the  noise  field  were  completely  uncorrelated,  these  curves  (see  figures  C1-C7) 
would  be  straight  lines  at  10  log  N,  where  N  is  the  number  of  sensors.  For  this  data  there 
were  5  sensors,  so  10  log  5  =  7  dB  for  uncorrelated  noise.  As  can  be  seen  by  examining  the 
sequence  of  figures  Cl  through  C7,  the  least  gain  is  realized  at  D/E  =  0  (figure  Cl),  and 
progressively  greater  array  gain  results  from  steering  beams  progressively  farther  up,  or 
down,  from  horizontal  (D/E  =  0).  The  term  D/E  indicates  the  depression/elevation  angle 
with  positive  angles  up  from  the  horizontal  and  negative,  angles  down  from  the  horizontal. 
Only  beginning  at  D/E  -  30  does  the  optimal  array  gain  curve  begin  to  exceed  the  10  log  N 
(7  dB)  level.  The  anomalous  peak  in  the  optimal  gain  curve  at  105  Hz  was  due  to  noise 
radiated  from  the  nearby  surface  vessel,  which  at  this  time  was  only  700  feet  from  the  sus 
pension  buoy.  Eor  the  eloser-iu  targets  which  are  expected  to  have  arrival  angles  up  to  60° 
D/E  (figure  (5).  improvements  of  up  to  15  dB  are  possible  using  adaptive  beamforming  at 
frequencies  above  80  llz.  Eor  conventional  beamforming,  improvements  of  only  3  or  4  dB 
less  are  possible. 

(U)  Fora  theoretical  description  of  array  gain  the  reader  is  referred  to  references  1 

and  2. 

(U)  Figures  D1  through  D42  present  calculated  array  gains  for  the  same  simulated 
noise  fields  as  for  the  directional  response  computations  above,  showing  significant  gains 
which  were  possible  even  from  this  array  without  close  attention  to  array  lilt,  vertical  mo¬ 
tion,  and  array  tether. 


CONFIDENTIAL 


13 


CONFIDENTIAL 


Frequency  Versus  Ik-pression/Flevation  Angle  Plots 

(U)  An  overall  view  of  the  recorded  noise  field  is  presented  in  the  I'Requency  vs 
Ik-pression/Flevation  angle  < l-'K I )/li)  plots.  These  plots  are  collections  of  single  line-bvarii  g 
responses  at  many  frequencies.  All  of  the  lines  on  a  plot  are  normalized  to  a  common  maxi¬ 
mum  power  point;  the  contour  lines  represent  areas  of  equal  power.  On  the  FRD/F  plots 
presented  in  this  report,  these  contours  are  1-dB  apart.  Figure  3  is  the  FRD/H  plot  for  the 
adaptive  processor.  The  dynamic  range  of  the  ploi  is  approximately  30  dll  with  the  maxi¬ 
mum  power  response  occurring  at  about  0=8 5  (near  horizontal)  and  frequency  =  76  Hz. 

The  noise  field  above  50  Hz  is  fairly  directional.  T  he  power  for  0  near  90°  appears  13  dB 
higher  than  the  powers  at  0-0°  and  0  -  1X0°.  The  harmonic  content  of  noise  at  0  =  85°  is 
also  revealed  by  the  adaptive  bearing  response.  (Note  that  the  dark  areas  near  0°  (straight 
up)  and  180°  (straight  down)  at  frequencies  above  150  Hz  are  dips  in  the  power  response.) 
The  apparent  isotropic  nature  of  the  noise  field  below  30  Hz  is  artificial;  it  is  caused  by  the 
bandpass  filters  applied  to  the  data  (as  mentioned  earlier).  Figure  4  is  the  FRD/H  plot  lor 
the  conventional  (timc-shift-und-siim)  processor.  Its  dynamic  range  is  approximately  16  dl> 
with  a  maximum  response  at  0  ~  72“  and  frequency  -*  76  Hz.  Almost  all  of  the  dark  areas 
on  this  plot  are  dips  rather  than  peaks.  Also,  the  frequency  content  present  in  the  conven¬ 
tional  FRD/H  is  not  as  well  defined.  The  field  does  not  appear  so  directional  on  the  conven¬ 
tional  FRD/K  as  on  the  adaptive.  There  is  only  a  10  dB  separation  between  power  at  0  =  90° 
and  power  at  0  ~  0°  and  0  =  180°.  This  generally  agrees  with  the  single  frequency  results 
(figures  A 1  through  A1 2). 

(U)  Figures  A 1  through  A 1 2  wore  generated  from  the  actual  data  at  frequencies 
selected  from  the  FRD/E  contour  plots  The  same  array  configuration  was  used  in  several 
simulated  noise  fields  at  approximately  the  same  frequencies.  Figures  Bl  through  B72  were 
generated  from  simulated,  vertically-direetional  noise  fields.  Comparison  of  the  actual  data 
figures  with  the  latter  figures  allows  conclusions  regarding  the  directionality  of  the  actual 
noise  field. 

(U)  Figures  Bl  through  B72  show  the  potential  resolving  power  as  calculated  from 
the  array  geometry  (at  the  frequencies  indicated)  by  plotting  bearing  vs  power  responses  for 
conventionally  and  adaptively  processed  beamforming  with  simulated,  directional  noise 
fields.  These  simulated  noise  fields  were  obtained  from  Naval  Air  Development  Center 
(NADC).  The  simulated  noise  field  is  the  solid  curve  on  the  graph,  while  the  output  of  Un¬ 
conventional  processor  is  the  curve  with  the  squares  and  the  output  of  the  adaptive  proces¬ 
sor  is  the  curve  with  the  circles.  -40  dB  of  uneorrelated  noise  was  added  to  each  of  the 
simulated  noise  fields.  As  stated  earlier,  the  unusually  high  level  of  uneorrelated  circuit 
noise  apparently  prevented  the  actual  array  implementation  from  realizing  its  potential  for 
resolving  directionality. 

(U)  The  analysis  output  which  quantifies  the  improvement  that  can  be  realized  by 
using  a  vertical  array  of  sensors  compared  with  using  only  a  single  omnidirectional  sensor  is 
ARRAY  (IAIN  vs  FREQUENCY.  Ail  example  of  this  output  is  shown  for  the  January  1970 
Mediterranean  data  (as  described  earlier)  as  figures  Cl  through  C7. 

(U)  Figures  1)1  through  1)42  show  the  array  gains  for  the  simulated  noise  fields. 

-40  dB  of  uneorrelated  noise  was  added  to  each  of  the  simulated  noise  fields. 
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13  abstract 


(U)  Data  were  taken  from  a  5-sensor  vertical  array  in  the  Ionian  Basin.  The  data  were  reduced  via  sophisticated 
analytical  processes  to  find  the  vertical  directionality  of  the  ambient  noise  field.  The  directionality  from  the  measured 
data  was  compared  to  several  noise  fields  simulated  a  priori.  The  measured  directionality  was  found  not  to  correspond 
well  to  any  of  the  a  priori  simulated  fields. 

(U)  It  was  observed  tnat  the  amount  of  uncorrelated  (circuit)  noise  lias  a  strong  influence  on  the  measured 
directionality.  Large  percentages  of  uncorrelated  noise  depress  the  ability  of  an  adaptively  processed  array  to  resolve  the 
field.  A  new  simulated  noise  field  was  generated  to  include  the  same  percentage  of  uncorrclatcd  noise  as  was  reported 
for  the  recording  instrumentation  (!0  percent).  This  simulation  was  found  to  compare  very  closely  with  the  actual  array 
response.  Since  the  high  uncorrclatcd  noise  probably  is  a  result  of  the  recording  rather  than  the  actual  field,  the  direc¬ 
tionality  was  again  computed  after  10  percent  uncorrclated  noise  was  removed  from  the  spectral  cross  power  matrix 
form  ot  the  actual  data.  The  resultant  directionality  was  then  compared  with  that  from  a  simulated  field  will  -m  uncor- 
rclatcd  noise,  which  is  presented  as  a  candidate  directionality  model  for  the  Mediterranean  Sea  ambient  noise  at  97  Hz. 

(U)  This  report  is  presented  as  a  prototype  example  for  directionality  analysis  of  data  from  n.  .hi-cicmcnt 
arrays  (VLAM,  FLU’,  ACODAC,  MO  DANS,  SONODIVER,  1TASS,  etc.)  Adaptive  bcamformin0;  ABE)  algorithms  were 
employed  as  the  analysis  tools.  Array  gain  versus  frequency  results  arc  also  shown  for  actual  ,md  simulated  noise  fields. 
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